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Direct detection experiments for neutralino dark matter in the Milky Way 
are examined within the framework of SUGRA models with R-parity invari- 
ance and grand unification at the GUT scale, Mq- Models of this type ap- 
ply to a large number of phenomena, and all existing bounds on the SUSY 
parameter space due to current experimental constraints are included. For 
models with universal soft breaking at Mq (mSUGRA), the Higgs mass and 
b — > 57 constraints imply that the gaugino mass, my 2, obeys my 2 >(300- 
400) GeV putting most of the parameter space in the co-annihilation domain 
where there is a relatively narrow band in the mo — rny 2 plane. For fi > 
we find that the neutralino -proton cross section ~ 10~ 10 pb for m\i% < 1 
TeV, making almost all of this parameter space accessible to future planned 
detectors. For fi < 0, however, there will be large regions of parameter space 
with cross sections < 10~ 12 pb, and hence unaccessible experimentally. If, 
however, the muon magnetic moment anomaly is confirmed, then \x > and 
mi/2 ~ 800 GeV. Models with non- universal soft breaking in the third gener- 
ation and Higgs sector can allow for new effects arising from additional early 
universe annihilation through the Z-channel pole. Here cross sections that will 
be accessible in the near future to the next generation of detectors can arise, 
and can even rise to the large values implied by the DAMA data. Thus dark 
matter detectors have the possibility of studying the the post-GUT physics 
that control the patterns of soft breaking. 
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1 Introduction 



The recent BOOMERanG, Maxima and DASI data has allowed a relatively precise de- 
termination of the mean amount of dark matter in the universe, and these results are 
consistent with other astronomical observations. Within the Milky Way itself, the amount 
of dark matter is estimated to be 

p DM = (0.3 - 0.5)GeV/cm 3 (1) 

Super symmetry with R-parity invar iance possesses a natural candidate for cold dark mat- 
ter (CDM), the lightest neutralino, Xii an d SUGRA models predict a relic density con- 
sistent with the astronomical observations of dark matter. Several methods for detecting 
the Milky Way neutralinos exist: 

(1) Annihilation of Xi m t ne halo of the Galaxy leading to anti-proton or positron 
signals. There have been several interesting analyses of these possibilities |], , but there 
are still uncertainties as to astronomical backgrounds. 

(2) Annihilation of the x? i n the center of the Sun or Earth leading to neutrinos and 
detection of the energetic by neutrino telescopes (AMANDA, Ice Cube, ANTARES). 
Recent analyses [0, |J indicate that these detectors can be sensitive to such signals, but for 
the Minimal Supersymmetric Standard Model (MSSM) one requires m^o > 200 GeV (i.e. 
m 1 / 2 > 500 GeV) and tan/3 > 10, and for SUGRA models one is restricted to tan/3 > 35 

I- 

(3) Direct detection by scattering of incident xi on nuclear targets of terrestrial de- 
tectors. Current detectors are sensitive to such events for Xi ~ P cross sections in the 
range 

a^_ p ~ 1 x 10~ 6 pb (2) 

with a possible improvement by a factor of 10 - 100 in the near future. Future detectors 
(GENIUS, Cryoarray, ZEPLIN IV) may be sensitive down to (10~ 9 — 10~ 10 ) pb and we 
will see that this would be sufficient to cover the parameter space of most SUGRA models. 

In the following we will consider SUGRA models with R-parity invariance based on 
grand unification at the GUT scale Mq = 2 x 10 16 GeV. In particular, we will consider 
two classes of models: Minimal supergravity models (mSUGRA || ^]) with universal soft 
breaking masses at Mq, and non- universal models with non universal soft breaking at Mq 
for the Higgs bosons and the third generation of squarks and sleptons . Here the gaugino 
masses (mi/2) and the cubic soft breaking masses (Aq) at Mq are assumed universal. 

SUGRA models apply to a wide range of phenomena, and data from different ex- 
periments interact with each other to greatly sharpen the predictions. We list here the 
important experimental constraints: 

Higgs mass: >114 GeV 0. The theoretical calculation of still has an an error of 
~ 3 GeV, and so we will (conservatively) interpret this bound to mean (theory) > 111 
GeV. 
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b — > S7 branching ratio. We take a 2a range around the central CLEO value |§: 

1.8 x 10~ 4 < B(B -> X s7 ) < 4.5 x 10~ 4 (3) 
Xi relic density: We assume here 

0.02 < n DM h 2 < 0.25 (4) 

The lower bound takes into account of the possibility that there is more than one species 
of DM. However, results are insensitive to raising it to 0.05 or 0.10. 
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Figure 1: Corridors in the mo — vnxii plane allowed by the relic density constraints for 
(bottom to top) tan/3 = 10, 30, 40, A = and fi > 0. The lower bound on mu2 is due to 
the mfc lower bound for tan/3 = 10, due to the b — > 37 bound for tan/3 = 40, while both 
these contribute equally for tan/3 = 30. The short lines cutting the channels represent 
upper bound from the g^ — 2 experiment. [17] 

Muon a M = (g tl — 2)/2 anomaly. The Brookhaven E821 experiment || reported a 2.6a 
deviation from the Standard Model value in their measurement of the muon magnetic 
moment. Recently a sign error in the theoretical calculation |L0j, 111] has reduced this to 



a 1.6a anomaly, though recent measurements |L2] used to calculate the hadronic contri- 
bution may have raised the deviation. Since there is a great deal of more data currently 
being analyzed (with results due this spring) that will reduce the errors by a factor of 



^2.5, we will assume here that there is a deviation in a M due to SUGRA of amount 

llxlO- 10 <af GRA <75xKr 10 



(5) 



We will, however, state our results with and without including this anomaly. 

To illustrate how the different experimental constraints affect the SUSY parameter 
space, we consider the mSUGRA example: 
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Figure 2: Corridors in the mo — mi« plane allowed by the relic density constraint for 
tan/3 = 40, fi > and (bottom to top) A = 0, — 2m 1 / 2 , 4my 2 . the curves terminate at 
the lower end due to the b — > 37 constraint except forAo = 4mi/2 which terminates due 
to the rrih constraint. The short lines cutting the corridors represent the upper bound on 
mi/2 due to the g M — 2 experiment. [17] 



(1) The rrih and b —>■ sj constraints put a lower bound on mi/2' 

mi/2 ~ (300 - 400)GeV (6) 

which means m^o ~ (120 — 160) GeV (since m^o = 0.4m!/ 2 ). (2) Eq.(6) now means 
that most of the parameter space is in the fi — x? co-annihilation domain in the relic 
density calculation. Then mo (the squark and slepton soft breaking mass) is approximately 
determined by rrii/2 as can be seen in Figs. 1 and 2. (3) If we include the a M anomaly, since 
g sugra - g a d ecreas i n g function of m^j and mo , the lower bound of Eq.(5) produces an 
upper bound on m^j and the positive sign of a M implies that the /1 parameter is positive. 
In addition one gets a lower bound on tanbeta of tan [3 > 5. Thus the parameter space 
has begun to be strongly constrained, allowing for more precise predictions. In order to 
carry out detailed calculations, however, it is necessary to include a number of analyses 
to obtain accurate results. We list some of these here: 

Two loop gauge and one loop Yukawa renormalization group equations (RGE) are used 
in going from M G to the electroweak weak scale M EW , and QCD RGE are used below 
Mew f° r the light quark contributions. Two loop and pole mass corrections are included 



in the calculation of m^. One loop corrections to m?, and m T [13|, [14]] are included which 
are important at large tan f3. Large tan (3 NLO SUSY corrections to b — >• 37 [IB] are 
included. In calculating the relic density, all stau-neutralino co- annihilation channels are 
included, and this calculation is done in a fashion valid for both small and large tan j3. 
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Figure 3: &x°-p for mSUGRA for // < 0, A = 1500 GeV, for tan/3 = 6 (short dash), 
tan/5 = 8 (dotted), tan/? = 10 (solid), tan/3 = 20 (dot-dash), tan/? = 25 (dashed). Note 
that the tan/3 = 6 curve terminates at low rai/2 due to the Higgs mass constraint, and 
the other curves terminate at low m 1( / 2 due to the b — > 57 constraint [18]. 

We do not include Yukawa unification or proton decay constraints, since these depend 
sensitively on post-GUT physics, about which little is known. 

2 mSUGRA MODEL 

The mSUGRA model is the simplest, and hence most predictive of the supergravity models 
in that it depends on only four new parameters and one sign (in addition to the usual 
SM parameters). We take these new parameters to be mo and mi/2 (the universal soft 
breaking scalar and gaugino masses at Mq), Aq (the universal cubic soft breaking mass 
at Mq) , tan/5 =< H 2 > / < Hi > at the electroweak scale (where < H 2 > gives rise to 
up quark masses and < Hi > to down quark masses) and the sign of \x (the Higgs mixing 
parameter which appears in the superpotential as fiHiH 2 ). We examine these parameters 
over the range mo,mi/ 2 < 1 TeV, 2 < tan/? < 50, \A \ < 4mi/ 2 . The bound on rrii/ 2 
corresponds to the gluino mass bound of m- g < 2.5 GeV which is also the reach of the 
LHC. 

The relic density analysis involves calculating the annihilation cross section for neu- 
tralinos in the early universe. This characteristically proceeds through Z and Higgs s- 
channel poles (Z, h, H, A where H and A are heavy CP even and CP odd Higgs bosons) 
and through t-channel sfermion poles. However, if there is a second particle which be- 
comes nearly degenerate with the neutralino, one must include it in the early universe 
annihilation processes, which then leads to the co-annihilation phenomena. In mSUGRA 
models, this accidental near degeneracy occurs naturally for the light stau, fi. One can 
understand this semi-quantitatively by considering the low and intermediate tan /3 region 
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Figure 4: o"x?-p as a f unc tion of the neutralino mass m^o for tan j3 = 40, fi > for 
Aq = — 2m 1 / 2 , Ami/2, from bottom to top. The curves terminate at small m^o due to the 
b — > S7 constraint for v4 = and -2mi/2 and due to the Higgs mass bound (m^ > 114 
GeV) for A = 4mi/2. The curves terminate at large rn^o due to the lower bound on a M 
ofEq. (5) [17]. 



where the RGE give for the right selectron, e^, and the neutralino the following masses 
at the electroweak scale: 

m \ R = ml + 0.15m 2 /2 - sin 2 6 w cos 2 (3 (7) 
m|o = 0.16m 2 /2 (8) 

the numerics coming from the RGE analysis. The last term in Eq. (7) = (40 GeV) 2 . 
Thus for m = the cr will become degenerate with the x? a t ^1/2 = 400 GeV, and 
co-annihilation thus begins at mi/ 2 = (350 — 400) GeV. As mi/ 2 increases, tuq must be 
raised in lock step (to keep mg H > m^o). More precisely, it is the light stau, which is the 
lightest slepton that dominates the co-annihilation phenomena. However, one ends up 
with corridors in the mo — mi/2 plane for allowed relic density with m closely correlated 
with mi/2 increasing as mi/ 2 does, as seen in Figs. 1 and 2. 

For dark matter detectors with heavy nuclei targets, the spin independent neutralino - 
nucleus cross section dominates, which allows one to extract the Xi -proton cross section, 
cr^o_p. The basic quark diagrams for this scattering go through s-channel squark poles 
and t-channel Higgs (h, H) poles. The general features of cr^o_ p that explain its properties 
are the following: 

cr^o_ p increases with increasing tanbeta (9) 

0^0 _ p decreases with increasing mi/2 and increasing mo (10) 

Since co-annihilation generally correlates m and m^, if mi/2 increases so does m (at 
fixed tanbeta and A ). 
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The smallest cross sections occur for the case /x < 0. This is because a special can- 
cellation can occur over a fairly wide range of tan/3 and [19|, [18] driving the cross 
section below 1CT 13 pb. This is illustrated in Fig. 3. In these regions, there would be no 
hope for currently planned dark matter detectors to be able to detect Milky Way neu- 
tralinos. However, if the a M anomaly is confirmed by the new BNL E821 data (currently 
being analyzed), then /j < is forbidden, and the special cancelations do not occur for 
/i > 0. Large cross sections can then occur for large tanbeta. Thus is seen in Fig. 4 for 
tan/? = 40, with rrih > 114 GeV. If the Higgs mass bound were to rise, the lower bounds 
on mi/2 would increase. Thus for rrih > 120 GeV, one has m x o > (200, 215, 246) GeV for 
A = (-2, 0, 4)mi /2 . 

The lowest cross sections for ii > are expected to occur for small tanbeta and large 
mi/2- This is seen in Fig. 5 for tan/? = 10 where one also sees that decreasing A$ gives 
smaller cross sections. In general one finds 



o'xl-p ~ 10 10 pb for fi > 0, mi/2 < 1 TeV 
Such cross sections are within the reach of future planned detectors. 
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Figure 5: a x°-p as a f unc tion of m x o for tan/3 = 10, /i > 0, > 114 GeV for Aq = 
(upper curve), Aq = — 4mi/2 (lower curve). The termination at low m^o is due to the rrih 
bound for A = 0, and the b — > 57 bound for A = —4mi/ 2 
is due to the lower bound on a M of Eq. (5) [17]. 
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3 NON-UNIVERSAL MODELS 

New results can occur if we relax the universality of the squark, slepton and soft breaking 
Higgs masses at Mq- To maintain the flavor changing neutral current bounds, we do this 
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Figure 6: Effect of a nonuniversal Higgs soft breaking mass enhancing the Z° s-channel 
pole contribution in the early universe annihilation, for the case of 62 =1, tan/5 = 40, 
A = mi/2, > 0. The lower band is the usual f\ coannihilation region. The upper 
band is an additional region satisfying the relic density constraint arising from increased 
annihilation via the Z° pole due to the decrease in /1 2 increasing the higgsino content of 
the neutralino[18]. 



only in the third generation and for the Higgs bosons, 
breaking masses at Mq as follows: 



One may parameterize the soft 



ml(l + 6 1 ] 
[l + Ss, 



ITLr 



mJl + S 6 



711 



m 



m 



Hz 

2 

2 



rrir 



(I + 62 
[1 + 5, 



m 



TR 



rrir 



m 2 Jl + 5 7 



(12) 



with — 1 < 5i < +1. While the non- universal models introduce a number of new parame- 
ters, it is possible to understand qualitatively what effects they produce on dark matter 
detection rates, since the parameter /z 2 governs much of the physics. Thus as u 2 decreases 
(increases), the higgsino content of the neutralino increases (decreases), and then cr^o_ p 
increases (decreases). One can further see semi-quanitatively the dependence of u 2 on the 
non-universal parameters for low and intermediate tan/3 where the RGE may be solved 



analytically [20 
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(13) 



where t = tan/? and Dq = 1 — (m t /200sin/3) 2 . In general Dq is small i.e. Dq = 0.25, and 
one sees that the universal part of the m 2 , contribution is quite small, and it does not take 
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Figure 7: 



a^o _ p as a function of mi/2 



0.4mi/ 2 ) for tan/3 = 40, // > 0, m h > 114 



GeV, A = 7711/2 for 5 2 = 1. The lower curve is for the fj — %5 co-annihilation channel, 
and the dashed band is for the Z s-channel annihilation allowed by non-universal soft 
breaking. The curves terminate at low rai/2 due to the b —>■ 37 constraint. The vertical 
lines show the termination at high 777.1/2 due to the lower bound on of Eq. (5) [21]. 



a great deal of non-universal contribution to produce additional effects. 

Most interesting things happen when /i 2 is decreased, since the increased Higgsino 
content of the neutralino increases the Xi~Xi~Z coupling, and this coupling opens a 
new annihilation channel through the Z-pole in the relic density calculations. As a simple 
example we consider the case where only the H 2 soft breaking mass is affected i. e. 5 2 = 1 
and all other 5i = 0. Fig. 6 shows the new allowed region in the 777.0 — 7774/2 plane for 
tan p = 40, A = 777,1/2 , /i > 0, and Fig. 7 shows the corresponding effect on the neutralino 
- proton cross section. On sees that the co-annihilation corridor is significantly raised and 
widened due to the new Z-channel annihilation, and the cross section is significantly 
increased. The next round of upgraded dark matter detectors should be able to reach 
parts of this parameter space if such a non-universality were to occur. 

As a second example we consider a soft breaking pattern consistent with an SU(5) 
invariant model with Siq(= 63 = 64 = 85) = —0.7, and all other 5, = 0. Here the tr soft 
breaking mass is reduced , i.e m| = 777,q(1 + 5 5 ) < m^. Thus the f\ — Xi co- annihilation 
occurs at a larger value of m than in mSUGRA. In addition again a new Z-channel 
neutralino annihilation channel occurs since /1 2 is reduced. The effects are shown in Figs. 
8 and 9 for tan/3 = 40, A = 777.1/2, yU > 0. Again the cross sections are larger, and should 
be accessible to CDMS when it moves to the Soudan mine and to GENIUS. 

The maximum value of cr^_ p for fixed tan/3 and A Q occurs when we chose the non- 
universalities to minimize /z 2 . This occurs when #1,3,4 < and 62 > 0. This is shown 
in Fig. 10 where the maximum cross section is plotted for A = 0, tan/3 = 12 (upper 
curve), tan/3 = 7 (lower curve). The bound that rrih > 114 GeV, eliminates the region 
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Figure 8: Allowed regions in the m — 7^1/2 plane for the case tan/5 = 40, A = m^, 
\i > 0. The bottom curve is the mSUGRA fx coannihilation band of Fig. 1 (shown for 
reference). The middle band is the actual f\ coannihilation band when 5 10 = —0.7. The 
top band is an additional allowed region due to the enhancement of the Z° s-channel 
annihilation arising from the nonuniversality lowering the value of fi 2 and hence raising 
the higgsino content of the neutralino. For m 1( / 2 ~ 500 GeV, the two bands overlap [18]. 

with m^o < 100 GeV. However, one sees for this case that it is possible to have detection 
cross sections in the region of the DAMA data. 

4 CONCLUSIONS 

We have discussed here direct detection of Milky Way neutralinos for SUGRA type models 
with R-parity invariance and grand unification at the GUT scale. By combining data 
from a variety of sources, e.g. Higgs mass bound, b — > 57 branching ratio, relic density 
constraints and the possible new muon magnetic moment anomaly of the BNL E821 
experiment, one can greatly sharpen predictions. 

For the mSUGRA model, the rrih and b — *■ S7 bounds create a lower bound on m 1( / 2 
of m 1/2 ~ (300 - 400) GeV (i. e. m^o ~ (120 - 140) GeV). Thus puts the parameter 
space mostly in the 77 — Xi co-annihilation domain, which strongly correlates m with 
mi/2- For fi > and myi < ITeV, one finds cr^o_ p ~ 10 -10 pb which is within the upper 
reach of future planned dark matter detectors, while for /j, < there will be large regions 
unaccessible to such detectors. If the a M anomaly is confirmed, then \i > and < 800 
GeV. 

Non-universal soft breaking models allow one to raise o"^o_ p by a factor as large as 
10 - 100, which could account for the large cross sections of the DAMA data. They can 
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Figure 9: 0"x5-p as a function °f m i/2 for tan/3 = 40, /i > 0, A = m-i/2 and m fc > 114 
GeV. The lower curve is for the bottom of the fi — x? co-annihilation corridor, and the 
upper curve is for the top of the Z channel band. The termination at low is due 
to the b — > 57 constraint, and the vertical lines are the upper bound on due to the 
lower bound of a M of Eq. (5) [21]. 



also open new allowed regions of the mo — mx/2 plane from the Z channel annihilation 
in the relic density calculation. The new Z-channel regions have larger cross sections, 
though still below the DAMA region, but they should be accessible when CDMS is in the 
SOUDAN mine and to the GENIUS-TF detector. Thus dark matter detectors should be 
able to investigate the nature of SUSY soft breaking, i.e. the nature of the post-GUT 
physics that determine the soft breaking pattern. 
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